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Heat  exchangers  and  pre-reformers  are  critical  devices  for  high  temperature  fuel  cell  systems.  It  is  rec¬ 
ommended  to  incorporate  a  compact  heat  exchange  and  the  pre-reformer  when  considering  the  limited 
space  and  cost.  The  volume-resistance  characteristic  modeling  technique  is  introduced  here  to  meet  the 
requirement  for  quick  dynamic  and  real  time  simulations.  The  distribution  characteristics  along  the  heat 
exchange  reformer  length  direction  are  presented,  and  some  key  effect  factors  are  studied.  The  tran¬ 
sient  behaviors  are  investigated  for  different  step-change  conditions,  such  as  mass  flow  rate  and  inlet 
temperature.  This  can  provide  some  references  and  tools  for  the  fuel  cell  system  design  and  optimization. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

High  temperature  fuel  cell  systems  are  an  attractive  emerg¬ 
ing  technology  for  stationary  power  generation,  especially  for  the 
distributed  generation  [1].  Today,  there  are  mainly  two  types  of 
high  temperature  fuel  cell  systems,  including  the  molten  carbonate 
fuel  cell  (MCFC)  and  solid  oxide  fuel  cell  (SOFC),  which  are  gener¬ 
ally  operated  at  high  temperatures  ranging  from  823  K  to  1273  K. 
Several  advantages  of  this  setup  are  listed  in  Ref.  [2].  The  main 
advantages  of  both  fuel  cells  are  related  to  what  could  be  done 
with  the  waste  heat  and  how  they  can  be  used  to  reform  fuels,  pro¬ 
vide  heat,  and  drive  engines.  Therefore,  high  temperature  fuel  cell 
systems  can  never  be  simply  considered  as  fuel  cells;  instead,  they 
must  always  be  thought  of  as  an  integral  part  of  a  complete  fuel 
processing  and  heat  generating  system  [2]. 

Steam  reforming  is  a  well-established  industrial  fuel  process  for 
producing  hydrogen  or  synthetic  gas  from  natural  gas,  other  hydro¬ 
carbon  fuels,  and  alcohols  [3].  In  the  high  temperature  fuel  cell 
systems,  the  pre-reformer  is  usually  needed  for  fuel  processing. 
Due  to  the  high  endothermic  reaction,  a  great  amount  of  heat  must 
be  provided  from  the  outside,  such  as  waste  heat  from  the  fuel  cell, 
catalyst  combustion,  etc. 
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High  temperature  heat  exchangers  are  widely  used  in  the  high 
temperature  fuel  cell/gas  turbine  system,  closed  cycle  gas  turbine 
system,  high  temperature  gas  cooled  reactors,  and  other  thermal 
power  systems.  It  is  an  effective  method  of  improving  the  whole 
system  efficiency  [4].  Compact  heat  exchangers  are  generally  char¬ 
acterized  by  extended  surfaces  with  large  surface  area/volume 
ratios  that  are  often  configured  in  either  plate-fin  or  tube-fin 
arrangements  [5].  In  a  plate-fin  exchanger,  many  augmented  sur¬ 
face  types  are  used:  plain-fins,  wavy  fins,  offset  strip  fins,  perforated 
fins,  pin  fins,  and  louvered  fins.  Offset  strip  fins,  which  have  a  high 
degree  of  surface  compactness  and  feasible  manufacturing,  are  very 
widely  applied. 

In  general,  the  high  temperature  heat  exchanger  is  used  to  pre¬ 
heat  the  air  or  fuel,  while  the  pre-reformer  is  used  to  produce 
hydrogen  rich  fuel  from  methane  or  other  hydrocarbons.  Fig.  1 
shows  one  of  the  fuel  cell  systems,  which  consists  of  a  direct 
internal  reforming  solid  oxide  fuel  cell  (DIR-SOFC),  a  high  temper¬ 
ature  heat  exchanger  (HTHE),  a  low  temperature  heat  exchanger 
(LTHE),  a  pre-reformer,  a  gas  turbine,  a  generator,  etc.  In  order  to 
simplify  the  system,  reduce  the  cost,  and  improve  the  fuel  cell 
system’s  efficiency,  it  is  suggested  that  a  compact  heat  exchange 
reformer  replace  the  heat  exchanger  and  the  pre-reformer.  The 
new  fuel  cell  system  is  illustrated  in  Fig.  2.  The  offset  strip  fin  heat 
exchanger  and  pre-reformer  are  combined  into  the  heat  exchange 
reformer.  In  this  device  with  the  counter-flow  type,  the  high  tem¬ 
perature  waste  gas  from  the  fuel  cell  flows  in  the  hot  passage, 
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Nomenclature 

A  area(m2) 

C  molar  concentration  (mol  m-3 ) 

Cp  specific  heat  capacity  (kj  kg-1  K-1 ) 

Dh  hydraulic  diameter  (m) 

DEN  parameter  used  in  Table  1 

/  fanning  friction  factor 

G  mass  flow  rate  (kg  s-1 ) 

Gm  mass  velocity  (kg  m-2  s-1 ) 

J  Colburn  factor 

k  parameters  used  in  Table  1 ,  or  geometry  parameter 

used  in  formula  (8)  (m) 

I<  parameter  used  in  Table  1 

l  offset  strip  fin  length  (m) 

L  heat  exchanger  length  (m) 

M  molecular  weight  (kg  mol-1 ) 

n  number 

p  partial  pressure  of  component  i  in  the  cold  fuel  pas¬ 

sage  (Pa) 

P  pressure  (Pa) 

Pr  Prandtl  number 

R  gas  constant  (J  mol-1  K-1 ) 

Re  Reynolds  number 

S  passage  heat  transfer  surface  (m2) 

St  Stanton  number 

t  fin  or  plate  thickness  (m),  time  (s) 

T  temperature  (K) 

u  velocity  (ms-1) 

U  wet  perimeter  (m) 

W  whole  heat  exchanger  width  (m) 

X  passage  width  (m) 

Y  passage  height  (m) 

Greek  letters 

a  convective  heat  transfer  coefficient  (kj  m-1  s-1  K-1 ) 

or  dimensionless  geometry  parameter  used  in  for¬ 
mula  (13)  and  (15) 

y  dimensionless  geometry  parameter  used  in  formula 

(13) and  (15) 

8  dimensionless  geometry  parameter  used  in  formula 

(13) and  (15) 

AH,  A H°  enthalpy  change  and  enthalpy  change  at  the  stan¬ 
dard  state  (kj  mol-1 ) 

A P  pressure  loss  (Pa) 

p  fin  efficiency 

X  thermal  conductivity  (kj  m-1  s-1  K-1 ) 

/i  dynamic  viscosity  (Pas) 

p  density  (kg  m-3) 

a  friction  resistance 

Subscripts 
c  cold  side 

f  fin 

h  hot  side 

i  fuel  component 

w  solid  fin  structure 

(I)  steam  reforming  reaction 

(II)  gas  shifting  reaction 

(III)  C02  direct  reforming  reaction 


and  the  fuel  flows  in  the  cold  passage.  In  particular,  the  Ni  cat¬ 
alyst  is  coated  on  the  fuel  passage  surface  [6,7].  When  the  fuel 
flows  along  the  passage,  the  endothermic  steam  reforming  reac¬ 
tion  will  take  place  using  the  heat  transferring  from  the  hot 
side. 

Several  kinds  of  compact  heat  exchange  reformers  have 
been  investigated  and  designed  by  previous  authors.  In  2001, 
Kawasaki  Heavy  Industries  in  Japan  developed  a  plate-fin  heat- 
exchange  reformer  with  highly  dispersed  catalyst  [8].  A  planar 
micro-channel  concept  was  proposed  by  Pacific  Northwest 
National  Laboratories  (PNNL),  but  this  kind  of  micro-channel 
device  is  oriented  toward  the  low  to  medium  power  range 
(20-500  W)  for  man-portable  applications  [9,10].  A  novel  micro 
fuel  processor  for  PEMFCs  with  heat  generation  by  catalytic 
combustion  was  developed  and  characterized  in  South  Korea 
[11-13]. 

All  these  previous  works  were  mainly  developed  based  on 
experiments,  but  the  steady  state  and  dynamic  performance 
simulations  have  not  been  investigated  in  detail.  The  heat  sup¬ 
plied  for  the  methane  steam  reforming  reaction  has  different 
sources,  such  as  catalytic  combustion  [11,12]  and  auto-thermal 
methane  reforming  reactions  [10].  The  purposes  are  mainly 
for  the  portable  devices  [9,10]  or  the  low  temperature  fuel 
cell  [11-13].  Here,  the  waste  heat  from  the  high  temperature 
fuel  cell  systems  will  be  used  as  the  heat  resource  in  the 
compact  heat  exchange  reformer  for  the  steam  reforming  reac¬ 
tion. 

This  paper  aims  to:  design  a  compact  heat  exchange  reformer 
for  the  high  temperature  fuel  cell  systems;  develop  a  real  time 
simulation  model  using  the  volume-resistance  characteristic  mod¬ 
eling  technique;  study  the  steady-state  distribution  characteristics 
by  considering  local  fluid  properties  such  as  pressure,  velocity, 
density,  heat  capacity,  thermal  conductivity,  dynamic  viscos¬ 
ity,  etc.;  discuss  some  factors  that  will  affect  the  performance 
of  the  reformer  during  steady-state  operation  under  the  same 
operating  condition;  and  finally,  investigate  dynamic  behavior 
under  different  input  parameters  including  step-change  condi¬ 
tions. 

2.  Mathematical  models 

2.1  Configuration 

The  configuration  of  the  heat  exchange  reformer  is  similar  to  the 
compact  heat  exchanger.  The  only  difference  is  that  the  catalyst  is 
coated  in  the  cold  passage  to  make  steam-reforming  reactions  take 
place. 

As  shown  in  Fig.  3,  the  configuration  of  the  offset  strip  fin  heat 
exchanger  is  adopted  here.  The  fin  surface  is  broken  into  a  number 
of  smaller  sections.  Generally,  each  type  of  fin  is  characterized  by 
its  width  X,  height  Y,  thickness  t,  and  length  of  the  offset  strip  fin 
l.  The  detailed  configuration  can  also  be  found  in  other  references 
for  the  heat  exchanger  [14-18]. 

Taking  the  hot  passage  as  an  example,  the  calculations  for  indi¬ 
vidual  geometry  variables  are  listed  as  following: 

W 

Passage  number  :  nh  =  - -  (1) 

Xh  +  th 

Offset  strip  number  :  nh/  =  (2) 

*h 

Cross-area  of  flow  passage  :  Ah  =  nhXhYh  (3) 


Heat  transfer  surface  of  flow  passage  : 


Sh  —  2tth(Xh  +  Yh)L  4-  Hh/Hh(Xh  +  Yh  +  th)th 


(4) 
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Fig.  1.  Schematic  view  of  the  traditional  SOFC/GT  hybrid  system. 


Wet  perimeter :  Uh  =  (5) 

Hydraulic  diameter :  Dh  =  (6) 

uh 

2.2.  Passage  fin  efficiency 

The  passage  fin  efficiency  rj0  is  given  by  Rosehnow  et  al.  [18]  as 

rio  =  1  -  j(l  -  i/f)  (7) 

where  the  secondary  heat  transfer  area  of  a  stream  5f  for  the  hot 
passage  equals  Sh.  The  total  area  of  the  heat  exchanger  S  is  cal¬ 


culated  by  the  sum  of  the  primary  heat  transfer  surface  and  the 
secondary  heat  transfer  area  of  a  stream. 

According  to  Rosehnow  et  al.  [15,18],  the  fin  efficiency  for  the 
offset  strip  fin  with  a  rectangular  section  can  be  approximated  by 

tanh(mhkh) 

mhkh  <8) 

where  mh  =  y^h^hAh/h.  kh  =  Yh/( 2  -  th). 

Finally,  the  fin  efficiency  can  be  simplified  by 

>?0,h  =  1  -  y  (1  ~>?f,h)  (9) 

xh  +  Yh 

The  fin  efficiency  is  mainly  influenced  by  the  material,  configu¬ 
ration  of  the  fin,  and  the  heat  transfer  coefficient  between  the  fin 
and  the  flow. 


Fig.  2.  Schematic  view  of  the  SOFC/GT  hybrid  system  with  novel  concept  heat  exchange  reformer. 


H.  Zhang  et  al.  /  Journal  of  Power  Sources  183  (2008 )  282-294 


285 


Table  1 

Reaction  and  its  rate  in  the  heat  exchange  reformer  (Xu  and  Froment  [19]) 


(I) CH4  +  H20<^C0  +  3H2 

(II) C0  +  H20^C02  +  H2 

(III) CH4  +  2H20  C02  +  4H2 


k  (  Ph  Pc o  \  1 

(16) Rm  =  (pcH4PH20  -  Kei  )  Blip 

(17) R(»)  =  sir  (PcoPh2°  _  P-3£r)  bsp 

(18) J?(hi)  =  i  (pch4P^0J  -  T^1) 


Table  2 

Hot  passage  dynamic  mathematical  model 


dph 

d 1 


9(Ph»h2 


(23) Mass  conservation  equation  : 

(24) Momentum  conservation  equation  : 

(25) Energy  conservation  equation  :  ^  =  -uh  ^  -  |^^(Th  -  Tw) 


9(Phuh) 

dx 

9(Phuh)  _ 
dt  - 


.  9Ph 

dx 


Uh«h 

Ah 


Fig.  3.  Flow  (a)  and  fin  (b)  diagram  of  heat  exchange  reformer. 


Table  3 

Cold  passage  dynamic  mathematical  model 


(26)Mass  conservation  equation  :  ^ 


(27)Momentum  conservation  equation  : 

djpcllc2)  , 


Ucac 

Ac 


(28)Energy  conservation  equation  : 

,1  3 Tc  ScOfc^O.C  fT  T  \  i  1 

Uc~dZ  -  PcCpcacL^c-'wJ+^Yc 


3(pcUc) 

dt 


dTc  _ 
dt  - 

E 


(-AH)kRk 


ke{(I),  (II),  (III)} 


2.4.  Heat  transfer  coefficient 

Generally,  the  heat  transfer  coefficient  a  is  related  to  the  Colburn 
factor  [15,17,18]  and  is  expressed  as 

a  =JGmCpPr~2/3  (14) 


2.3.  Pressure  loss 


The  frictional  pressure  loss  across  an  offset  strip  fin  passage  and 
at  any  associated  entry,  exit,  and  turning  loss  [15],  can  be  expressed 
by 


AP  =  4f 


(10) 


where  Gm  =  pu. 

Here,  turning  losses  are  neglected,  so  the  pressure  loss  at  per 
unit  length  can  be  expressed  by 


A P 

~T 


(M 


(11) 


Let  the  friction  resistance  a  =  1/2 (fpu2), 
then 

dP  _Ua 
d x~  A  ' 


(12) 


The  fanning  friction  factor  /  has  been  developed  by  many 
authors.  The  correction  equation  has  been  developed  as  the  fol¬ 
lowing  equation  by  Manglik  and  Bergles  [17]  from  the  data  of  Kays 
and  London  [14] 


where  the  Colburn  factor  J=StPr2l 3  and  the  Prandtl  number 
Pr=  fiCpIX. 

The  correction  equation  has  also  been  developed  by  Manglik  and 
Bergles  [17]  from  the  data  of  Kays  and  London  [14].  The  Colburn 
factor  expression  is  listed  as 

J  =  0.6522f?e-°-5403a(-0154150149V“00678 

x[l  +  5.269  X  lo-SRe1340^-504^0-456^-1  055]0'1  (15) 


2.5.  Steam  reforming 

In  the  cold  fuel  passage,  the  steam  reforming  reaction  (I), 
water  gas  shift  reaction  (II),  and  C02  direct  reforming  reactions  of 
methane  (III)  are  carried  out  over  a  Ni  catalyst  coat  on  the  passage 
surface  at  sufficiently  high  temperatures,  typically  above  773  K. 

Kinetic  rate  equations  for  the  reactions  (I)— (III)  are  adopted  from 
Xu  and  Froment  [19].  The  three  kinetic  rate  equations  are  listed  in 
Table  1  as  well.  The  enthalpy  changes  of  chemical  reactions  are 
calculated  according  to  the  author  Smith  et  al.  [20] 

AH(I)  =  AH(°I}  -  16373.61 


/  =  g tQ243Re~°  7422 a~0A856 8°  3053 y~°  2659 

x[l  +7.669  X  10-8Re4-429a0-92053.767K0.236]01 


+R  (  7.951  Tc  -  4.354e  -  3 T2  +  0.7213e  -  6 Tc3 


0.097e5  \ 
Tc  ) 
(19) 


(13) 
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Table  4 

Heat  exchange  reformer  volume-resistance  characteristic  model 


Hot  passage 

(30)^  = 


L*h,l  1*11,1-^11,2 


/oi'idGh,2  A  ph,l-ph,2  IT  _ 

(31 )  nr  =  Ah 


(32) 


d  t 
dTh 


^*h,2  *h,l  ~  *h,2 


^h®h,  2  /"t’  t1  ^ 

dx  A^C^V^h,2  V21 


Cold  passage 

^  dQ,j,2  _  _  uc,2Q,i,2~uc,lQ,i,l 

dx 


(33  )T 

y  vi>kRk,2£  is{CH4,  H2,  CO,  C02,  H20| 


*«{<I),  (II),  (III)l 

fOA\dPc,2  _  RTc.2  CC2~Gc 
rit  ~  ] 


dt  Mc  Ac  dx 

(35)%1=4cM1-Uc<tc,1 


=  ** 


Gc,i  Tc,\~TC,2  SCac,i  (T 

pC;1CpclAcL^c,l 


Tw’l)+  Pc,lCpcl  Yc 


(-AH)k>1Uk>1 


dx 

cPc,i  'c 

ke{(I),(H),  (HI)} 

Solid  structure 

f'in\ dlw  2  v  Tw  3-2Tw  2+Tw  i 

(3?  )-ar  ='f* - ^ - 

®h,2^h  Cl -  y  \  ®C,2^C  /y  y  \ 

MwCnw(7h>2  JVfwCn  (ic’2  iw>2^ 


Table  6 

Key  simulation  parameters  under  the  basic  condition 


Simulation  conditions 

Cold  fuel 

Inlet  mass  flow  rate  (kg  s-1 ) 

0.06 

Inlet  temperature  (I<) 

898 

Fluid  molar  fraction 

0.25  CH4,  0.75  H20  (STC  =  3:1) 

Outlet  pressure  (Pa) 

1.0E  +  5 

Hot  waste  gas 

Inlet  mass  flow  rate  (kg  s-1 ) 

0.4 

Inlet  temperature  (K) 

1200 

Fluid  molar  fraction 

0.1  C02,  0.2  H20,  0.1  02,  0.6  N2 

Outlet  pressure  (Pa) 

1.0E  +  5 

2.6.  Mathematical  model 

To  simplify  the  complexity  of  the  mathematical  model,  some 
assumptions  [4,21  ]  adopted  in  the  theoretic  analysis  are  presented 
as  follows: 

(1)  The  heat  exchange  reformer  is  adiabatic  to  the  surrounding; 

(2)  The  viscosity  dissipation  effects  are  neglected; 

(3)  The  parameters  are  considered  to  be  uniform  over  a  cross- 
section,  one  dimensional  flow  along  the  passage,  without  inside 
circumfluence; 


A H(II)  =  AH°n)  -  7756.56  +  R  (l  .86TC  -  0.27e  -  3 Tc2  +  1,1^4g5) 

(20) 

AH(|||)  =  AH(°m)  -  26125.07 

+R  (l0.6577c-4.624e-3Tc2  +  0.7213e-6Tc3  +  1~°^7e5^ 

C  (21) 


Table  5 

Geometry  and  properties  parameters  of  heat  exchange  reformer 


System  geometry  parameters 

Length 

1  m 

Width 

0.5  m 

Height 

0.532  m 

Hot  passage 

Width 

4.5E-3  m 

Height 

6.5E-3  m 

Offset  strip  fin  length 

0.05  m 

Fin  thickness 

3.0E-3  m 

Cold  passage 

Width 

4.5E-3  m 

Height 

5.0E-3  m 

Offset  strip  fin  length 

0.05  m 

Fin  thickness 

5.0E-3  m 

Separator 

Thickness 

1.0E-3  m 

Solid  structure  properties  (SiC  ceramic  [27-29]) 

Density 

3100  kg  irr3 

Heat  capacity 

0.640  kj  kg-1  K-1 

Thermal  conductivity 

0.080  kjm-1  s-1  K-1 

Catalyst  properties 

Thickness 

5.0E-5  m 

Density 

2355  kg  irr3 

Catalyst  reduced  activity 

0.003 

Dimensionless  properties  Dimensionless  properties  Dimensionless  pressure 
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(4)  For  the  horizontal  fluid,  the  effect  of  height  change  can  be  omit¬ 
ted. 

In  the  cold  fuel  passage,  the  chemical  species  are  CH4,  H2,  CO, 
C02,  and  H20.  Species  mass  balances  in  the  cold  fuel  passage  are 
considered 

3Cc,i  _  ,  3Cc,i 
at  c  dx 

+  Vi’kRki  ie(CH4’  H2 ’  CO,  C02,  H20} 

ke{(I),  (II),  (III)) 

(22) 


The  mass,  momentum,  and  energy  conservation  equations  for 
the  hot  passage  and  cold  passage  are  established  in  Tables  2  and  3, 
respectively.  In  the  hot  passage,  the  heat  transfer  to  the  solid 
structure  is  considered.  Due  to  the  very  thin  catalyst  coat,  the 
enthalpy  changes  of  the  reactions  (I)— (III)  are  also  considered  in 
the  cold  passage,  in  addition  to  the  heat  transferred  from  the  solid 
structure. 

For  the  solid  structures,  such  as  the  fins  and  the  separators,  the 
temperature  is  considered  to  be  uniform  at  the  same  cross-section. 
The  energy  conservation  equation  is  written  as 


3Tw  „32Tw  , 

~dF=,<^  + 


^h^h*?0,h 

MwCPw 


(Tw- 


Th)  + 


ucScmic 

MwCPw 


(Tw- 


Tc) 


(29) 


Fig.  9.  STC  effect  on  the  methane  (a)  and  hydrogen  (b)  molar  fraction  distributions. 
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Fig.  11.  CRA  effect  on  the  methane  (a)  and  hydrogen  (b)  molar  fraction  distributions. 


The  heat  conductivity  coefficient  is  I<=LXwAwIMwCpw, 
the  cross-area  of  solid  structure  is  Aw  =  2Wt+nh(Xh  +  Yh  + 
th)th  +  nc(Xc  +  Yc  +  tc)tc,  and  the  mass  is  Mw  =  PwAwL. 

The  control  equations  of  the  heat  exchange  reformer  are  strongly 
coupled.  In  addition  to  the  partial  differential  equations  presented 
above,  two  perfect  state  equations  P=f[pJ)  for  the  hot  and  cold 
passages  are  also  needed  in  order  to  compose  a  close  equation 
set. 

3.  Simulation  modeling  and  conditions 

3.1.  Volume -resistance  characteristic  model 

In  general,  nonlinear  partial  differential  equations  will  be  solved 
using  a  numerical  difference  method,  which  makes  it  very  dif¬ 
ficult  to  obtain  theoretical  solutions.  This  kind  of  equation  set 
is  very  simple  and  clear.  However,  stability  is  one  crucial  fac¬ 
tor  in  using  a  difference  algorithm.  In  addition,  the  time  step  for 
the  difference  algorithm  is  usually  very  short,  so  it  is  very  time 
consuming  [4]. 

In  order  to  avoid  the  coupled  iteration  between  the  flow  rate 
and  pressure,  the  volume-resistance  characteristic  modeling  tech¬ 
nique  [4,22]  is  introduced  into  the  heat  exchange  reformer.  This 
modeling  technique  is  based  on  the  lumped-distributed  parameter 
method,  which  can  obtain  a  set  of  ordinary  differential  equations 
from  partial  differential  equations. 


The  volume-resistance  characteristic  model  is  listed  in  Table  4 
in  detail. 

3.2.  Simulation  conditions 

In  addition  to  the  configuration  and  geometry  parameters 
of  the  heat  exchange  reformer,  as  shown  in  Table  5,  and  fluid 
properties  calculated  at  the  local  position,  some  boundary  con¬ 
ditions  were  also  required  to  carry  out  the  simulation.  These 
included  inlet  flow  rate,  fluid  composition,  and  the  inlet  tem¬ 
perature  and  outlet  pressure  of  both  the  hot  and  cold  streams 
(Table  6). 

At  the  same  time,  some  conditions  are  used  to  solve  the  equa¬ 
tions;  for  example,  the  thermal  flux  density  of  both  the  solid 
structure  inlet  and  outlet  are  considered  to  be  zero.  As  a  result, 
contrasted  to  the  center  difference  algorithm  in  the  middle  of  the 
solid  structure,  the  difference  algorithms  for  both  the  front  and  end 
modules  are  treated  independently. 

4.  Results  and  discussion 

In  this  section,  due  to  the  high  cost  of  the  complicated  exper¬ 
iments,  only  simulation  studies  are  employed  on  a  counter-flow 
type  heat  exchange  reformer.  Section  4.1  provides  the  dis¬ 
tributed  characteristics  of  some  important  parameters,  such  as 
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fuel  species,  temperature,  and  fluid  properties  (pressure,  den¬ 
sity,  velocity,  heat  capacity,  thermal  conductivity,  and  dynamic 
viscosity),  under  steady-state  conditions.  Section  4.2  compares 
and  analyzes  the  results  under  different  input  parameter  con¬ 
ditions,  such  as  steam  to  carbon  ratio,  catalyst  reduced  activity, 
and  operating  outlet  pressure.  In  Section  4.3,  the  dynamic 
behaviors  of  the  compact  heat  exchanger  reformer  are  investi¬ 
gated. 

4.1.  Steady-state  result  analysis 

For  the  rated  condition,  some  related  parameters  are  presented 
in  Table  6,  such  as  inlet  temperature,  mass  flow  rate,  molar  fraction, 
and  outlet  pressure. 

Fig.  4  presents  the  fuel  molar  fraction  along  the  heat  exchange 
reformer  length.  The  flow  direction  in  the  fuel  channel  is  from 
1.0  to  0  in  the  figures,  so  all  the  parameters  in  the  fuel  chan¬ 
nel  should  be  understood  to  proceed  from  1.0  to  0.  At  the  cold 
fuel  passage  inlet,  the  fluid  only  contains  methane  and  water.  The 
steam  reforming  reaction  takes  place  on  the  surface  of  the  catalyst 
along  the  flow  direction.  Therefore,  the  methane  is  gradually  con¬ 
sumed.  The  methane  and  water  concentration  decreases  along  the 
flow  direction.  The  concentration  of  produced  hydrogen  gradually 
increases.  The  methane  steam  reforming  reaction  has  two  simulta¬ 
neous  effects.  The  carbon  monoxide  molar  fraction  increases  and 


Fig.  13.  Cold  passage  outlet  pressure  effect  on  the  methane  (a)  and  hydrogen  (b) 
molar  fraction  distributions. 


the  carbon  dioxide  molar  fraction  increases  along  the  flow  direc¬ 
tion.  At  the  exit,  the  flow  composition  is  4.24%  of  CFI4, 45.35%  of  FI2, 
10.00%  of  CO,  3.84%  of  C02,  and  36.57%  of  H20. 

The  temperature  profiles  of  the  cold  stream,  hot  stream,  and 
solid  structure  along  the  heat  exchange  reformer  length  are 
presented  in  Fig.  5.  Because  of  the  high  endothermic  methane 
reforming  reaction,  the  cold  fuel  temperature  decreases  a  lit¬ 
tle  at  the  entrance.  Then,  the  cold  fuel  temperature  increases 
along  its  flow  direction  due  to  the  heat  transfer  from  hot 
gas.  The  temperatures  of  the  hot  gas  stream  and  the  solid 
structure  decrease  along  the  heat  exchange  reformer  length.  It 
should  be  noted  that  the  temperature  curve  is  just  the  line 
between  measured  points,  so  it  cannot  indicate  the  trend  at  both 
ends. 

The  pressure  profiles  in  the  cold  fuel  and  hot  gas  passages  are 
illustrated  in  Fig.  6.  Owing  to  the  friction  of  the  passage,  the  pressure 
loss  is  about  0.08%  in  the  cold  fuel  passage,  and  about  4.23%  in 
the  hot  gas  passage.  The  primary  reason  that  the  pressure  loss  is 
greater  in  the  hot  gas  passage  is  that  the  mass  flow  rate  in  the  hot 
gas  passage  is  larger  than  that  in  the  cold  passage.  Of  course,  the 
geometrical  configuration  is  a  key  factor  as  well. 

The  dimensionless  fluid  properties  (such  as  density,  velocity, 
heat  capacity,  thermal  conductivity,  and  dynamic  viscosity)  of  the 
cold  fuel  and  hot  gas  along  the  heat  exchange  reformer  are  illus- 
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Fig.  14.  Cold  passage  outlet  pressure  effect  on  the  cold  fuel  (a)  and  hot  gas  (b) 
temperature  distributions. 


trated  in  Figs.  7  and  8,  respectively.  The  dimensionless  properties 
are  defined  as  the  ratio  of  local  values  and  corresponding  inlet 
values,  which  can  be  calculated  by  the  inlet  conditions  in  the 
methods  depicted  in  Ref.  [23].  Examples  of  this  include  situations 
where:  the  density  is  based  on  the  gas  state  equation;  the  veloc¬ 
ity  is  calculated  by  the  mass  flow  rate,  density  and  the  channel 
cross-area;  the  heat  capacity  of  the  multi-component  gas  mix¬ 
ture  is  related  to  the  single  component  heat  capacity  and  the 
corresponding  molar  fraction;  the  dynamic  viscosity  of  the  multi- 
component  gas  mixture  is  based  on  the  Reichenberg’s  expression; 
the  thermal  conductivity  of  multi-component  gas  mixtures  is  based 
on  Wassiljewa’s  expression  and  the  Mason  and  Saxena  modifica¬ 
tion. 

The  density  is  related  to  the  pressure  and  the  temperature, 
which  are  decided  by  the  gas  state  equation  P=pRT.  In  the  cold 
fuel  passage,  the  temperature  increases  and  the  pressure  decreases, 
so  the  density  decreases  along  the  flow  direction  while,  in  the  hot 
gas  passage,  both  the  pressure  and  the  temperature  decrease.  The 
ratio  of  pressure  and  temperature  along  the  passage  is  increased, 
so  the  density  of  the  hot  gas  increases  along  the  flow  direc¬ 
tion. 

Two  primary  factors  that  affect  the  velocity  are  the  mass 
flow  rate  and  the  density.  Here,  the  mass  flow  rate  is  con¬ 
stant,  and  the  velocity  is  mainly  determined  by  the  density.  That 
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Fig.  15.  Dynamic  response  of  the  temperatures  at  the  cold  (a)  and  hot  (b)  passage 
exits  when  cold  fuel  mass  flow  rate  up  by  10%. 


is  to  say,  the  velocity  increases  in  the  cold  fuel  passage  and 
decreases  in  the  hot  gas  passage,  following  the  trend  of  the  den¬ 
sity. 

Specific  heat  capacity,  thermal  conductivity,  and  dynamic 
viscosity  are  primarily  influenced  by  the  temperature  and 
the  gas  composition.  This  has  been  discussed  by  Todd  and 
Young  [24]  and  Wang  et  al.  [22]  for  high  temperature 
SOFCs. 

4.2.  Analysis  of  the  influence  of  some  parameters 

In  this  section,  some  key  parameters  that  affect  the  heat 
exchange  reformer  performance  are  investigated,  such  as  the  steam 
to  carbon  ratio  (STC),  catalyst  reduced  activity  (CRA),  and  passage 
operating  pressure. 

4.2 A.  Steam  to  carbon  ratio 

In  general,  the  STC  must  be  greater  than  2.0  to  avoid  carbon 
coking  in  the  fuel  lines,  reformer,  and  fuel  cell  stack  [25].  The  effect 
of  different  STCs  on  the  heat  exchange  reformer  is  presented  in 
Figs.  9  and  10. 

Fig.  9  presents  effect  of  STC  on  the  methane  and  hydrogen  distri¬ 
bution  along  the  heat  exchange  reformer.  In  the  internal  reforming 
high  temperature  fuel  cell,  the  endothermic  reforming  reaction  will 
cause  a  great  temperature  gradient,  which  could  decrease  the  life 
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Fig.  16.  Dynamic  response  of  methane  (a),  hydrogen  (b)  and  water  (c)  distributions  when  cold  fuel  mass  flow  rate  up  by  10%. 


of  the  fuel  cell  stack  due  to  excessive  thermal  stress.  Therefore, 
too  much  remaining  methane  would  be  no  good  for  the  steady 
operation  of  the  high  temperature  fuel  cell.  With  the  STC  chang¬ 
ing  from  2:1  to  4:1,  less  methane  remains  at  the  exit  (Fig.  9(a)), 
while  the  hydrogen  molar  fraction  at  the  exit  is  almost  the  same  as 
at  the  entrance  (Fig.  9(b)).  Therefore,  a  suitable  and  acceptable  STC 
is  essential  for  the  internal  reformation  of  high  temperature  fuel 
cells. 

The  temperature  distribution  of  cold  fuel  and  hot  gas  is  illus¬ 
trated  in  Fig.  10.  When  the  STC  changes  from  2:1  to  4:1,  less 
methane  is  provided  at  the  inlet,  and  less  heat  is  needed  for  the 
steam  reforming  reaction.  Meanwhile,  a  higher  STC  will  result 
in  a  higher  rate  of  the  exothermic  water  gas-shift  reaction,  so 


the  temperature  curves  of  both  the  cold  and  hot  stream  are 
higher. 

4.2.2.  Catalyst  reduced  activity 

The  CRA  is  defined  as  the  ratio  between  the  activity  of  the  cata¬ 
lyst  in  use  and  that  of  a  conventional  Ni  catalyst  (Xu  and  Froment 
[19])  at  typical  feed  conditions  (temperature,  pressure,  and  compo¬ 
sition)  [26].  The  CRA  is  the  key  factor  in  determining  the  reforming 
reaction  rate.  For  the  rated  case,  the  CRA  is  defined  as  0.003  [7] 
in  Table  5.  Figs.  11  and  12  present  the  effect  of  the  CRA  on  the 
performance  of  the  heat  exchange  reformer. 

The  influence  on  the  methane  and  hydrogen  molar  fraction  dis¬ 
tribution  along  the  heat  exchange  reformer  is  shown  in  Fig.  11. 


Cold  passage  exit  temperature-K 
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Fig.  17.  Dynamic  response  of  the  temperatures  at  the  cold  (a)  and 
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passage  (b)  exits  when  the  hot  inlet  temperature  down  to  1100 1<. 


Fig.  18.  Dynamic  response  of  methane  (a),  hydrogen  (b)  and  water  (c)  distributions  when  the  hot  inlet  temperature  down  to  1100 1<. 


H.  Zhang  et  al.  /  Journal  of  Power  Sources  183  (2008 )  282-294 


293 


When  the  CRA  changes  from  0.0015  to  0.006,  the  rate  of  the 
methane  reforming  reaction  increases,  so  more  methane  is  con¬ 
sumed  (Fig.  11(a))  and  more  hydrogen  is  produced  (Fig.  11(b)).  More 
heat  is  needed  to  satisfy  the  requirements  of  the  high  endothermic 
reaction,  so  the  temperature  curves  of  both  the  cold  and  hot  stream 
are  lower  (Fig.  12). 

4.2.3.  Passage  operating  pressure 

The  passage  pressure  often  changes  with  the  operation 
condition,  even  during  malfunctions  or  damage.  The  effect 
of  the  cold  passage  outlet  pressure  on  the  heat  exchange 
reformer  is  investigated  in  this  section  and  illustrated  in 
Figs.  13  and  14. 

The  cold  passage  outlet  pressure  has  little  influence  on  the  heat 
exchange  reformer  performance.  When  the  passage  pressure  is  ele¬ 
vated  from  lE+5Pa  to  4E+5Pa,  less  methane  is  consumed,  less 
hydrogen  is  produced  (Fig.  13),  and  less  heat  is  needed  for  the 
methane  steam  reforming  reaction,  so  the  cold  fuel  and  hot  gas 
temperatures  are  higher  (Fig.  14). 

4.3.  Dynamic  simulation  result 

In  this  section,  the  transient  behaviors  of  the  compact  heat 
exchange  reformer  are  investigated.  Several  step-change  input 
parameters  (such  as  inlet  mass  flow  rate  and  inlet  temperature  of 
both  the  cold  and  hot  stream)  are  imposed  when  the  device  has 
been  operated  for  500  s. 

Fig.  15  illustrates  the  dynamic  response  of  the  temperatures  at 
the  cold  and  hot  passage  exits,  when  the  cold  fuel  mass  flow  rate 
has  a  step  increase  of  10%.  The  cold  passage  exit  temperature  has 
a  sudden  decrease  at  the  initial  period  due  to  the  step  input.  Then, 
because  of  the  great  thermal  inertia  of  the  solid  structure,  the  tem¬ 
perature  decreases  gradually.  Therefore,  the  temperature  at  the  cold 
passage  exit  decreases.  Owing  to  a  greater  cold  fuel  mass  flow  rate, 
more  heat  is  provided  from  the  hot  side,  so  the  temperature  at  the 
hot  passage  exit  has  a  gradual  decrease. 

Fig.  16  shows  the  dynamic  effect  on  methane,  hydrogen,  and  the 
water  molar  fraction  distribution  when  the  cold  fuel  mass  flow  rate 
has  a  step  increase  of  10%.  The  methane  and  water  molar  fraction 
increase  a  little,  while  the  hydrogen  decreases  a  little.  It  can  be 
shown  that  the  molar  fraction  has  a  little  change  when  the  cold 
fuel  inlet  mass  flow  rate  changes. 

Fig.  17  presents  the  dynamic  response  of  the  cold  fuel  and  hot 
gas  temperatures  when  the  hot  gas  inlet  temperature  decreases 
to  1100K  from  1200  K.  The  temperature  at  the  cold  passage  exit 
is  influenced  by  the  thermal  capacity  of  the  solid  structure,  and 
decreases  gradually.  Owing  to  the  decrease  of  the  inlet  tempera¬ 
ture,  the  temperature  at  the  hot  gas  passage  exit  also  undergoes 
a  decrease  (Fig.  17(b)).  When  the  temperature  of  the  cold  stream 
decreases,  the  rate  of  the  steam  reforming  reaction  will  be  slower. 
Therefore,  less  fuel  is  reformed,  which  can  be  shown  from  the 
methane  molar  fraction  distribution  in  Fig.  18(a);  less  hydrogen  is 
produced  (Fig.  18(b))  and  more  water  remains  (Fig.  18(c)). 

Based  on  all  the  dynamic  performance  figures  from  Figs.  15-18, 
the  inertial  delay  time  of  this  kind  of  heat  exchange  reformer  is 
about  3000  s.  Such  a  substantial  thermal  inertia  can  seriously  influ¬ 
ence  the  whole  fuel  cell  hybrid  system  transient  performance  and 
the  design  of  the  control  system. 

5.  Conclusions 

A  compact  heat  exchange  reformer  for  high  temperature  fuel  cell 
systems  is  presented  in  this  paper.  Based  on  the  volume-resistance 
characteristic  modeling  technique,  the  distributed-lumped  param¬ 
eter  method,  and  the  modular  modeling  idea,  a  simulation  model 


that  is  suited  for  quick  and  real  time  simulations  is  completed.  The 
model  can  predict  the  key  distribution  characteristic  parameters 
and  the  influence  of  some  factors,  such  as  the  steam  to  carbon 
ratio,  catalyst  reduced  activity,  and  passage  pressure.  The  dynamic 
results  indicate  that  this  kind  of  heat  exchange  reformer  has  a  great 
thermal  inertia. 

Both  the  model  and  modeling  method  will  be  useful  and  valu¬ 
able  for  other  heat  exchange  reformer  designs  and  optimization;  it 
can  also  provide  a  reference  for  the  design  of  the  control  system  in 
the  future. 
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